Some remarks on the visible points of a lattice 
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Abstract 

We comment on the set of visible points of a lattice and its Fourier trans- 
form, thus continuing and generalizing previous work by Schroeder Q and 
Mosseri 0. A closed formula in terms of Dirichlet series is obtained for the 
Bragg part of the Fourier transform. We compare this calculation with the 
outcome of an optical Fourier transform of the visible points of the 2D square 
lattice. 

Recently, Mosseri has given a nice and elegant description of the set of visible 
points of a lattice (i.e., those points except the origin that connect to the origin via 
a straight line without hitting any other lattice point in between) and some of its 
properties ||, continuing previous work by Schroeder [|l], [|. Both authors are most 
interested in the Fourier transform of this set. Unfortunately, the results of these 
two attempts are contradictory. The algebraic approach of [§] is essentially correct, 
but contains a couple of mistakes and unnecessary restrictions. In what follows, we 
add several remarks to straighten this out and then compare the 2D case with an 
optical experiment. 

Let A = Z&i © ■ ■ ■ © Z6 n be a lattice in nD space with linearly independent basis 
vectors b\, . . . , b n . The set F\ of visible points 0, [3|, H can be characterized as 

F\ = {mih H h m n b n \ gcd{m x , . . . , m n ) = 1}, (1) 

where gcd denotes the greatest common divisor. This set does not include the origin 
with respect to which it is defined. We follow || for the notation as far as possible. 

The set F\ is non-periodic and is left invariant by the group of lattice automor- 
phisms, Aut(A), which is isomorphic to Gl(n, Z), the group of integer nxn matrices 
with determinant ±1. This is seen from M 6 Aut(A) transforming fundamental 
cells of the lattice to other fundamental cells and hence visible points to visible 
points j|, ||. As a consequence, the set of visible points admits precisely the same 
point symmetry as the lattice A itself, the transformations are not restricted to pure 
rotations. This can clearly be seen from Fig. 1 which shows the case of the square 
lattice. 



How frequent are visible points? If p denotes the probability of a lattice point to 
be visible (defined through a volume limit which exists, compare chapter 3.8 of 0), 
we also have (feN = {1,2,3,...}) the probabilities 



Fix e £-F A ) = 



(2) 



because £-F\ = {mi&i + ■ • ■ + m n b n \ gcd(mi, . . . , m n ) = £}. On the other hand, we 
have 



A = {0} U \j£-F A 



(3) 



which is the correct version of eq. (1) in ||. Note that (Q) is the union of pairwise 
disjoint sets[] which will be essential in what follows. 

The single point 0, however, does not matter in the calculation of the probabilities 
wherefore we get 

OO OO 

l = ^GfF A )= P ^-. (4) 
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We can thus write p by means of Riemann's (^-function as 
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where V denotes the set of primes @, Obviously, p = for n = 1 (precisely 
two points are visible here), and p rapidly approaches 1 with increasing n, compare 
Fig. 2. For even n > 2, the ^-function is known to be transcendental f| |§ (and 



(2^ 



2(2m)! 



Bi m I , where Bo m are the Bernoulli numbers, 



thus irrational) through <^(2m) 
compare ||. Also, ((3) is irrational |Tt)|, while this question is open for odd n > 3. 



Note that the nonperiodicity of F\ follows from the irrationality of p, but not vice 
versa. 

A different approach uses the 3rd Mobius inversion formula 0, f|] where one starts 
from the characteristic function hr of a set T, which is one for points of the set and 
zero otherwise, and rewrites Eq. (0) as 



h\ = /i{o} + hi.p A . 



(6) 



Convergence is no problem because, for any given point, the right hand side is a 
finite sum. Now, using the inversion formula for such series M one obtains 



h 



e-A\{o} 



(7) 



with the Mobius function fi(£), cf. [§]. Note the differences to eqs. (1-3) in which 
result from the little mistake in the description of the lattice A. It is important 
though because | Y^t=\^{^)\ does not converge: it diverges faster than N 1 ^ 2 0, [11 



1 This property would be destroyed by the (unconventional) addition of the origin to the set Fa- 
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This indicates that the calculation of the Fourier transform requires some care: 
it is not clear that the set of visible points is strictly almost periodic in the sense 
that the Fourier transform need not be a sum over ^-distributions (or Bragg peaks) 
only. In fact, for n = 1, this diffractive part vanishes due to the existence of only 2 
visible points and the Fourier transform is 2 cos(kx) and thus continuous. This limit 
must nevertheless be covered by the correct treatment. In JJ], this is hidden by an 
infinite normalization factor - the formula would give vanishing Fourier transform 
for n — 1. 

Let us now consider the so-called structure factor in more detail. The structure 
factor Sr of a discrete set T is the Fourier transform of 5-scatterers of equal strength 
on all points of T, i.e., 

S r (k) = [Yl 6 ( x ' ~ x ) e ' lkX ' dx ' = Yl e ~ lkX - ( 8 ) 

J x£Y x£Y 

If T has inversion symmetry (as A and F\ do), we can also write Sr{k) = Y2 x &r e%kx = 
Szer cos (kx). Now, with Eq. (^), we find 

oo 

= £ j>(*)fc.A\{0}(*) " e ikx , (9) 

x£A 1=1 

to be understood in the distribution sense. Nevertheless, due to lack of convergence, 
an interchange of the sums is a subtle business. Though a complete treatment is 
desirable here, it is far beyond the scope of this comment and therefore left to 
a forthcoming publication [JT^ . But in the spirit of ||, one can indeed extract the 
point-like contributions, henceforth called Bragg peaks (though we are talking about 
amplitudes rather than intensities here) because other contributions are arbitrarily 
small in comparison to the 5-peaks in the thermodynamic limit. 

Before we continue, let us remark that the structure factor Sp A (k) is periodic in 
k w.r.t. the reciprocal lattice A*: we can alternatively write 



J x&K 



(10) 



with the characteristic function hp A defined above. But then, Sp A {k) is the convolu- 
tion of S\(k) - which is periodic ! - with the Fourier transform of the characteristic 
function hp A . Consequently, Sp A (k) itself is periodic. Furthermore, since the set Fa 
of visible points is invariant under Aut(A), the corresponding property applies to 
its Fourier transform, which is thus invariant under Aut(27rA*) - the latter again 
being isomorphic with the group Gl{n,'L). These two arguments apply to the full 
structure factor (and also to n — 1), not only to its Bragg peaks. 

To describe the latter, we now consider the reciprocal lattice 27rA* with A* = 
{y | x-y 6 Z Vx G A} = 7hb\ + • • • + Z6*, b*-bj = 5ij, a Bragg peak sits on every 
k = 2n(ribl + • • • + r n 6*) with rational coefficients, i.e., Tj G Q. If we write = Pi/qi 
with coprime integers p^, qi (where we take qi > for convenience), the peaks have 
strength H = H(q ll ...,q n ), 

H(q±, . . . , q n ) = V j—^- , a = lcm(q 1 , ...,q n ), (11) 
L — ' (may 1 

m=l v ' 
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relative to the Bragg peaks of the lattice A. Here, 1cm denotes the least common 
multiple. In particular, qi = . . . = q n = 1 gives 



tf(W) = V^ = -i-, (12) 



which is again the frequency of visible points, while H is zero unless all arguments 
are square free, because the Mobius function vanishes for numbers that are divisible 
by a square. This explains the "dark" lines in the diagram of intensities (i.e., lines 
without spots - hence "white" in Fig. 3 while really dark in Fig. 4) which are given 
- in kinematic approximation - by the absolute squares of the amplitudes ([TT]) . 

Eq. (|Tl~D gives a closed formula for the amplitudes of the Bragg peaks in form of 
an infinite sum, but can still be simplified considerably. If we consider 

for £, a G N, we can use standard techniques from the treatment of Dirichlet series, 
compare ||, to evaluate the sum. If a — pi ■ . . . -p r is the product of pairwise distinct 
primes, one can easily show by induction that 

'<»■»■■■■•»> = i^SrV (14) 

and hence, with the multiplicativity of the Mobius function and Eq. (|T 



Inserting this into (|TTD finally gives, with a = lcm(qi, . . . , q n ), the amplitude function 

''iufe' (16) 

This way, the calculation of the amplitudes relative to the central amplitude is 
reduced to the evaluation of a finite product, compare Fig. 3 for the 2D example 
of the square lattice. There, a Bragg peak is represented by a dot with radius 
proportional to the amplitude (i.e., the area is proportional to the intensity). 

The Fourier transform is periodic, although the set itself obviously is not - it 
is not even quasiperiodic because the Fourier module (i.e., the Z-module generated 
by the positions of the 5-peaks) is only countably but not finitely generated. The 
key structure is thus contained in one fundamental domain of the reciprocal lattice. 
This is indeed clearly seen in Figs. 3 and 4. The latter was obtained by optical 
Fourier transform of a finite portion of F %2 on an optical bench. Here, we used some 
9000 points of radius 15 % of the lattice constant. They were prepared as a negative 
slide made from 24 x 36 mm 2 documentary film. The light was extracted from a 
laser beam, widened by a high quality lens to a parallel beam covering the entire 
slide. Although the image is not perfect (one should not expect rapid convergence 
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in view of the remarks after ([])) and still shows some standard optical errors, it is 
nevertheless amazing how clear the structure with the subtle invariance properties 
is recovered. It is correctly described by the approach via the Mobius transform 
while previous numerical attempts [|I|] seem to fail for reasons we were not able to 
unravel. 

Of course, one can now deal with generalizations like the set of points visible 
to two or more observers which resembles the above situation in many respects, 
although the relative position of the observers does play a significant role. Unfortu- 
nately, we could not find similarly nice closed expressions for the Fourier transform 
wherefore we drop further details here. 
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Figure 1: Some visible points of the 2D square lattice. 
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Figure 2: Relative frequency of visible points in a lattice as a function of the dimen- 
sion. 
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Figure 3: Bragg part of the Fourier transform of the square lattice according to (jig). 
Here, the radius of the dots representing the Bragg peaks are proportional to the 
amplitude. 
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Figure 4: Optical Fourier transform of some 9000 visible points of the square lattice 
Z 2 , see text for details. 
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This figure "diffrac.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/math-ph/9903046vl 



